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ABSTRACT 

If a gamma ray burst with strong UV emission occurs in a molecular cloud, there 
will be observable consequences resulting from excitation of the surrounding H2. The 
UV pulse from the GRB will pump H2 into vibrationally-excited levels which produce 
strong absorption at wavelengths A < 1650 A. As a result, both the prompt flash 
and later afterglow will exhibit strong absorption shortward of 1650 A, with specific 
spectroscopic features. Such a cutoff in the emission from GRB 980329 may already have 
been observed by Fruchter et al.; if so, GRB 980329 was at redshift 3.0 < z < 4.4. BVRI 
photometry of GRB 990510 could also be explained by H2 absorption if GRB 990510 is 
at redshift 1.6 < z < 2.3. The fluorescence accompanying the UV pumping of the 
H-2 will result in UV emission from the GRB which can extend over days or months, 
depending on parameters of the ambient medium and beaming of the GRB flash. The 
7.5-13.6 eV fluorescent luminosity is ~ 10 41 ' 7 ergs s _1 for standard estimates of the 
parameters of the GRB and the ambient medium. Spectroscopy can distinguish this 
fluorescent emission from other possible sources of transient optical emission, such as a 
supernova. 

Subject headings: galaxies: ISM - gamma rays: bursts - molecular processes - ISM: 
clouds - ISM: molecules 



1. Introduction 



At least some gamma ray bursts (GRBs) are accompanied by intense optical emission, as 
demonstrated by detection of a 9th magnitude optical transient coinciding with GRB 990123 (Ak- 
erlof et al. 1999) and optical afterglows associated with other GRBs (e.g., GRB 990510: Stanek et 
al. 1999; Israel et al. 1999). Since GRBs may be associated with star- forming regions (Paczyhski 
1998), it is of interest to consider observable phenomena which would be indicative of molecular 
gas in the vicinity of the GRB. Here we show that there may be both observable absorption and 
emission due to H2 within several parsecs of a GRB. 

The energy radiated in the optical-UV flash and afterglow can be substantial, and can have 
dramatic effects on interstellar gas and dust in the vicinity of the GRB. The hv > 13.6 eV emission 
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will photoionize the gas, and the pulse of optical radiation will vaporize dust grains out to substantial 
distances from the GRB (Waxman, Draine & Phinney 1999). The gamma-rays and hard X-rays 
emitted by the GRB will contribute to ionization of the nearby gas, but the UV and soft X-rays have 
the dominant effect because of the much larger number of photons, and much larger photoionization 
cross sections. 

Ultraviolet radiation will destroy H2 in the vicinity of the GRB, but before destruction the 
typical H2 molecule will be excited to a vibrationally excited state by UV pumping. While cold 
molecular hydrogen can absorb only at wavelengths A < 1110 A, this vibrationally-excited H2 
can absorb strongly at wavelengths as long as 1650 A. The characteristic wavelength-dependent 
absorption by this vibrationally-excited H2 should be imprinted on the spectrum of the optical-UV 
flash when it reaches distant observers. Subsequent UV emission from the "afterglow" would also 
be subject to this characteristic absorption, due to vibrationally-excited H2 generated by the initial 
UV pulse plus additional UV from the afterglow. 

In addition to imposing a characteristic absorption spectrum on the optical transient, the 
optical-UV flash will generate fluorescent UV emission from surrounding H2. Because of light- 
travel time delays, the fluorescent emission will appear to extend over perhaps tens of days. The 
emission is strong enough that it could be detectable, and it would have a characteristic spectrum 
which would distinguish it from other possible sources of optical-UV emission. 

2. Photoexcitation of H2 

Prior to the GRB, the surrounding H2 is almost entirely in the first two or three rotational 
levels (J = 0,1,2) of the ground vibrational level (v = 0) of the ground electronic state (X 1 ^^). 
To photoexcite to the B X E+ and C 1 !!^ states requires hv > 11.2 and 12.3 eV, respectively. To 
photoionize H2 out of X 1 E+(v = 0, J = 0) requires hv > 15.4eV. The cold H2 itself is therefore 
essentially transparent to hv < 11.2 eV photons - the dominant absorption for hv < 11.2 eV is due 
to dust mixed with the H2. 

The first UV photons to arrive will most likely photoexcite the H2 via Lyman or Werner band 
transitions to the B or C electronic states, or photoionize the H2 to . H2 which is photoexcited 
to the B or C states will decay (in ~ 10~ 9 s) back to the ground electronic state, but typically 
to a vibrationally-excited level (e.g., v = 5) with the rotational quantum number changed by 
A J = 0, ±2. The lifetimes of the vibrationally-excited levels (~ 10 6 s) are long compared to the 
timescale for photoexcitation or photoionization [~ 10~ 3 s - see eq. (2,6,5)], so that depopulation of 
the rovibrationally-excited levels of H2 will be primarily by UV photoexcitation and photoionization. 

If the luminosity of the GRB pulse in 7.5-13.6 eV photons is 10 49 L4g ergs s _1 , and the spectrum 
L v oc v~ x l 2 as expected from simple models for emission from the forward shock (see Waxman, 
Draine, & Phinney 1999), then at a distance 10 19 i?ig cm from the GRB the specific energy density 
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u v is given by 

i/w 1/ = 4x 10" 14 x (1000 A/A) 1/2 ergs cm" 3 , (1) 

where 

X = 1.1 x 10 13 L49#r 9 2 ( 2 ) 

is the intensity at 1000A relative to the Habing (1968) estimate for the local interstellar radiation 
field. 

Stimulated emission in the UV transitions is negligible provided \ ^ 10 19 , so we consider 
photoexcitation out of level i of to vibration-rotation states of B 4 X+ and C^IlJ, followed 

by spontaneous decay to level j of X X S+. Let Tji be the rate for photopumping out of level i 
into level j, and let q\ p< ^ and be the rates for photodissociation and photoionization out of 
level i. Spontaneous decay via quadrupole transitions can be neglected if % S> 10 4 , and collisional 
deexcitation can be neglected if ran/ cm~ 3 <C \- If 



T = 



(3) 



then the H2 abundances p; L = 2n(H2(f«, Ji))/riR evolve according to 

j t Pi = Y. T W- ( 4 ) 
3 

We consider the 299 bound states of H2 with J < 29, and construct photopumping rates as 
described by Draine & Bertoldi (1996), using Lyman and Werner band oscillator strengths and 
dissociation probabilities from Abgrall et al ( 1993a, b) and Roueff (1993). We neglect attenuation 
of the radiation field due to absorption by H2, H, or dust. The rate of absorption of Lyman or 
Werner band photons for H2 in the v = 1, J = 3 level (for example), is 

Cf exc \v = 1, J = 3) = 3.34 x 10- 10 * s- 1 , (5) 

with a probability pd ~ 0.15 of photodissociation. 

If Lyman continuum radiation is present, then the H2 photoionization cross section from Yan, 
Sadeghour, & Dalgarno (1998) yields a photoionization rate 

= 1.96 x 10- 10 X s-\ (6) 

which we take to be independent of rovibrational excitation. We will also consider the case where 
all the Lyman continuum has been absorbed, so that = 0. 

For the photoionization rate of (6), the probability of returning to a bound rovibrational 
state following a photon absorption event is a = (1 — Pd)^^^ / {C^^ ^ + Cj P ^) ~ 0.54, the 
probability of ultimately being destroyed by photodissociation (rather than photoionization) is 
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PdCj PeXC ^ I '{Pd(j PeXC ^ + Cj P ^) ~ 0.20, and the mean number of fluorescent photons emitted per H2 
destroyed by radiation is 

1 C (PI) 

N& - = n V tod] (pi) * L36 ' ( 7 ) 

including transitions to the vibrational continuum of X 1 ^^. If, however, the illuminating spectrum 

is cut off at the Lyman limit, then = 0, and iVfl. ~ 6.7 - about 5 times as many fluorescent 
photons as for the case where Lyman continuum is present. 

A detailed study of the H2 photoexcitation would require modelling the time-dependent spec- 
trum of the radiation as it advances through the cloud. For a preliminary estimate of the rovibra- 
tional excitation we neglect changes in the spectrum of the UV pulse, and simply assume the H 2 to 
be irradiated by the spectrum (1), with and without Lyman continuum radiation. We assume the 
H2 to be "cold" at t = 0, with pj(0) = 0.5 for (v, J) = (0,0) and (0,1), and we follow the evolution 
of the level populations Pj(t) in the presence of the radiation field (1) until negligible H2 remains 
(%2pj < 10~ 6 ). We compute the time-integrated level populations Pj = pjdt. In Figure 1 we 
show the normalized level populations (fij = Pj/Ylj-^j versus level energy Ej. We see that even 
with Lyman continuum radiation, the time-averaged level populations have ~ 50% of the H2 in 
vibrationally-excited levels; if photoionization is suppressed, however, ~ 70% of the H2 is in levels 
v>l. 



3. Photoabsorption by UV-pumped H 2 

Waxman, Draine & Phinney (1999) discuss the destruction of dust by the GRB flash, and 
estimate the dust destruction radius to be Rd ~ 2.5 x 10 19 cm for a ~10 sec pulse with luminosity 
L49 ~ 0.55, and density «h < 10 5 cm~ 3 . Since we will estimate below that the radial extent of 
UV-pumped H2 will be considerably smaller than R^, it seems likely that dust destruction will take 
place so rapidly (in the leading edge of the GRB optical pulse) that we can neglect the effects of 
dust absorption and scattering when considering the excitation and destruction of the H2 . 

Photoexcitation by the GRB optical pulse produces rovibrationally excited H2 with permitted 
absorption lines in the 912 — 1650 A region. Figure 2 shows the absorption, smoothed to a resolution 
R = A/AA = 100, produced by column densities 7V(H 2 ) = 10 16 , 10 18 , and 10 20 cm" 2 with the 
distributions <p(v, J) from Figure 1. We assume Voigt profiles, approximated following Rodgers 
&z Williams (1974), with Doppler broadening parameter b = 3 km s^ 1 , for the 28765 permitted 
Lyman and Werner band lines between levels with J < 29. It is clear that the rovibrationally- 
excited H2 absorbs quite strongly for A < 1650 A (hi> > 7.5 eV): for iV(H2) ~ 10 18 cm~ 2 , it can 
be seen from Figure 2 that approximately 1/3 of the energy between 1650 and 912A has gone into 
photoexcitation of H2. In the discussion below we will assume that all of the photons between 1650 
and 912A are available for photoexcitation of H 2 . 
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The energy radiated in the afterglow may be comparable to or even exceed that emitted in 
a few tens of seconds around the emission peak. Let A^ on be the total number of hv > 13.6 eV 
photons, and N uv be the total number of 7.5-13.6 eV photons. For the F v ~ u^ 1 / 2 spectrum 
observed for the afterglow of GRB 990510 (Stanek et al. 1999) and expected from simple models 
(e.g., Waxman, Draine, & Phinney 1999), N uv /N ion = 0.35. 

Dust within ~ 2.5 x 10 19 cm of the GRB will be vaporized by hv < 7.5 eV photons in the 
initial flash, so dust absorption will be neglected for the hv > 7.5 eV photons of interest here. In 
an infinite cloud, the GRB will fully ionize the hydrogen and helium within a radius 

/ 771 \l/3 

R lon « 1.5 x 10 19 i-^-j cm , (8) 

where = nH/10 3 cm~ 3 , the total energy in 7.5-13.6 eV photons is 10 50 i?5o ergs, and we have 
assumed nHe/^H = 0.1. For r < i?j on we will suppose that 20% of the H2 is destroyed by photodis- 
sociation rather than by photoionization, (absorbing Afi. = 1.36 UV photons per photodissociation) 
thus absorbing a fraction 0.2N a X0.5N ioQ /1.2)/N uv = 0.33 of the 7.5-13.6 eV photons at r < R ion . 
The remaining 67% of the UV photons will photodissociate H2, with a probability pd ~ 0.15 per 
absorption, and therefore will dissociate the H2 out to a radius 

The column density of this photodissociation zone is 71h(-Rh 2 ~~ ^ion) = 4 x lO^ni^-E^g 3 cm" 2 . 

A more accurate treatment must simultaneously model photodissociation of the H2, photoion- 
ization of H, H2, and He, and attenuation of the radiation field by these absorption processes. 
This is particularly critical since we estimate Ru 2 ~ -Rion- For the present, however, we can safely 
conclude that not all the 7.5 — 13.6 eV photons which can pump H2 will have been absorbed by the 
radius R lon where all of the hv > 13.6 eV photons are exhausted. Therefore there will be a shell of 
rovibrationally excited H2 outside i?i on - Since the linear extent of the zone of partially-dissociated 
(and therefore rovibrationally excited) H2 is determined by the attenuation length for the 1650- 
912A photons which can photodissociate the H2, Fig. 2 would suggest that this layer of partially- 
dissociated (and rovibrationally-excited) H2 would have a column density N(H.2) ~ 10 20 cm -2 , and 
therefore a linear extent AR w2x lO 17 ^ 1 cm. Because the vibrationally-excited levels have life- 
times > 10 6 s against spontaneous decay, and collisional deexcitation is negligible for ns < 10, the 
vibrational excitation in the layer of partially-dissociated H2 will persist for weeks after the GRB. 

Attenuation by the rovibrationally-excited H2 will produce a strong "jump" near 1650 A in 
the spectrum of both the optical flash and the afterglow. The absorption is due to many distinct 
absorption lines, as is evident if the spectrum is shown with a resolution R = 1000, as in Figure 3. 

Fruchter (1999) found that the spectrum of the GRB 980329 optical transient showed a strong 
drop in flux between I (9000A) and R (6600A). Fruchter interpreted this drop as due to Lyman 
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a (A = 1215 A) absorption at redshift z ~ 5, with the Lyman forest then suppressing the flux at 
shorter wavelengths. However, recent deep imaging and spectroscopy of the host galaxy indicate 
that it must be at z < 5, and probably z < 4 (Djorgovski et al. 1999). If GRB 980329 occured in a 
molecular region, the observed drop in flux between I and R could have been due to the onset of H2 
absorption at 1650A in the host galaxy, in which case GRB 980329 was at redshift 3.0 < z < 4.4. 

Stanek et al. (1999) reported that VRI photometry of the afterglow of GRB 990510 were 
described by F u oc ^-0- 4 6±o.08^ but the B band flux falls below this prediction, indicating either a 
deviation from power-law behavior in the source, or additional extinction either near the GRB or 
in an intervening galaxy. Here we note that rovibrationally-excited H2 could account for this if the 
redshift of GRB 990510 is in the interval 1.6 < z < 2.3, so that the 1650 A absorption "edge" falls 
between B and V. This redshift for GRB 990510 is not inconsistent with the detection of absorption 
lines at z = 1.62 (Vreeswijk et al. 1999a,b). 



4. H2 Fluorescence 

As seen above, most of the 6.5 — 13.6 eV photons emitted by a GRB in a molecular cloud 
will be absorbed by H2. A fraction of the photoexcitations will result in photodissociation, but 
approximately 80% of the absorbed energy will be reradiated in UV fluorescent emission. 

Waxman, Draine & Phinney (1999) have shown that the prompt optical flash will result in 
destruction of the dust within ~ 10 pc of the GRB. If the observer is situated in the "beam" of the 
optical flash, then the fluorescent emission can reach the observer along a path from which the dust 
has been cleared, out to a distance Rd ~ 10 pc. There will be some attenuation of the fluorescent 
emission by H2 absorption, but such absorption is relatively weak for A > 1300 A (see Figure 2). 

If the prompt optical flash is radiated into a cone with half-width 9 ~ 0.1 (for a beaming 
factor 4tt/tt9 2 ss 400(0. I/O) 2 ), then the fluorescent emission from a shell of radius R will reach the 
observer spread out over a time Atj^fl. ~ R9 2 /c ~ 39(i?/10 19 cm)(#/0.1) 2 days. Because Atj^fl. 
is long compared to the GRB flash, for purposes of estimating the light curve of the fluorescent 
emission we can approximate the flash itself as a delta function. 

Suppose that a fraction f(r) of the H2 at radius r is destroyed by the prompt flash, with 
fluorescent energy e emitted per destroyed H2 molecule. If the observer is situated along the beam 
axis, then the fluorescent emission from regions with time delay < t is 

pct/0 2 poo 

E u (t)=ire u 9 2 f(r)n(R 2 )r 2 dr + TT€ct f(r)n(K 2 )rdr (10) 

Jo Jct/e 2 

If we now assume that n(H.2) is independent of radius, and f(r) = 1 for r < Rn 2 , and for r > Ru 2 , 
then the apparent luminosity is 

dE v TT 2 / C 2 t 2 \ 

(^)h 2 h. = = -^e v n(^ 2 )R^c 1 - -^-^ I (11) 
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We have computed the fluorescence spectrum ue u for H2 being destroyed by continuum radia- 
tion with I v oc i/ -1 / 2 , both with and without a cutoff at 912 A; the results are shown in Figure 4. 
The emitted radiation is concentrated in the 1650 — 912 A (7.5-13.6 eV) range. The values of ve v 
differ by a factor ~5, but according to our estimate, ~90% of the destroyed H is in the photoionized 
zone, so ve v ~ 1 x 10 -10 ergs/H2 is a representative average. There will be some "self-absorption" of 
this emission by the rovibrationally-excited H2 itself, an effect not included in the present analysis. 
If this fluorescent emission has also to traverse cold H2, it will suffer additional absorption, but 
only in a small number of absorption lines between 1110 and 912 A. 

The fluorescent emission estimated in eq. (12) appears to be a necessary consequence of a 
GRB in a molecular region. It is intriguing to note that the estimated luminosity and timescale 
are comparable to what would emerge from a supernova, although the spectrum of course is very 
dissimilar. 

Since the emission is concentrated between 1650 — 912 A, this process could not account for 
the emission plateau observed ~ 20 — 30 days after GRB 980326 (Bloom et al. 1999) unless the 
redshift z ss 3.5 (so that the 1600-1000 A emission would be redshifted to 7200-4500 A). However, 
at z = 3.5 the observed F u = OAfiJy at R would require vL v = 2 x 10 44 ergss -1 (for 17 = 0.3, 
Hq = 70km s _1 Mpc -1 ), or n l J 2 E^q 5000, much larger than anticipated for standard parameters 
(71,3 ss 1, E50 ps 1). Furthermore, the spectrum obtained by Bloom et al. 28 days after the GRB 
does not appear consistent with a redshifted version of Figure 4. For GRB 980326, then, it does 
not appear likely that the observed emission plateau is due to H2 fluorescence; the z ~ 1 supernova 
hypothesis proposed by Bloom et al. appears much more plausible. 

However, H2 fluorescence could perhaps be seen in other GRBs. GRBs at z ~ 1.5 — so that 
the 1600 A emission peak is redshifted to ~4000A - offer the best prospect for detection of this 
fluorescent emission, although it would still be quite faint for nominal parameters: (-^i/)max ~ 
2.5 x 10- 9 n 3 ( J R H2 /10 19 cm) 2 Jy at 4000A. 



5. Discussion 

The GRB pulse will create a large amount of H^ mixed with surviving H2 in the partially- 
dissociated shell. In these regions Hjj~ will be formed on a timescale ~ 10 6 (10 3 cm- 3 /n(H 2 )) s. This 
H3" will promptly react to produce other species, e.g., Kj" + CO — > HCO + + H2. The millimeter- 
wave lines of HCO + and other species might be observable in absorption against the afterglow, 
although this would probably require that the redshift first be determined by other methods (e.g. 
the H2 absorption lines discussed here). 
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Perna & Loeb (1998), Ghisellini et al. (1998), and Bottcher et al. (1998) have previously called 
attention to time- varying atomic and ionic absorption lines and X-ray absorption edges expected for 
GRBs, and Ghisellini et al. and Bottcher et al. pointed out that time-varying Fe K-a fluorescence 
might be observable from GRBs. Perna, Raymond & Loeb (1999) have pointed out that atomic 
and ionic emission lines might be observable from the photoionized remnants left behind by GRBs 
in relatively nearby galaxies. 

Observations of any of these phenomena would provide valuable clues to the nature of GRBs 
and their environments. If GRBs occur in molecular regions, the relatively strong H2 absorption 
spectrum may offer the best observational prospects for GRBs at redshift z > 1.3 so that A ss 1650 A 
is accessible to ground-based spectrographs. 

6. Summary 

The principal conclusions of this paper are: 

• Ultraviolet emission from a GRB will produce vibrationally-excited H2 which will produce 
strong line absorption for 1650 > A > 912 A. This absorption signature should be imprinted 
on the spectrum of the radiation reaching us from GRBs in molecular regions. 

• Absorption by vibrationally-excited H2 could be responsible for the pronounced drop in flux 
between R and I for GRB 980329 (Fruchter 1999), if GRB 980329 is at 3.0 < z < 4.4. 

• Absorption by vibrationally-excited H2 could account for the drop in flux between V and B 
for GRB 990510 (Stanek et al. 1999), if GRB 990510 is at 1.6 < z < 2.3 

• Ultraviolet emission from a GRB will produce fluorescent emission from whatever H2 may be 
present within tens of parsecs of the GRB. This fluorescence is potentially observable [see eq. 
(12) and Fig. 4]. 

I am grateful to Bohdan Paczyhski for many valuable discussions, including pointing out the 
possible application to GRB 980329, to S.G. Djorgovski for communicating results on GRB 980329 
in advance of publication, to E. Roueff for making available H2 wavelengths and oscillator strengths, 
and to Robert Lupton for availability of the SM software package. This work was supported in 
part by NSF grant AST-9619429. 
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Fig. 1. — Time- averaged fraction of the irradiated H2 in different vibration-rotation states, versus 
energy of the level, for an irradiating spectrum I u oc z^ 1 / 2 , with intensity 10 5 <C X ^ 10 19 . For 
a typical GRB flash with L49 ~ 1 we estimate x ~ 10 13 [ see ec L- (2)]- Collisional processes are 
neglected, valid for nu < 10~ 3 xcm -3 . Upper panel shows results for irradiation with cutoff at 
912 A; lower panel shows result if there is no break at the Lyman limit. 
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Fig. 2. — Transmission through a medium with H2 which has been rovibrationally excited by 
UV pumping by radiation with I u oc v~ x l 2 and an intensity 10 5 « x < 10 19 relative to the 
Habing intensity. Transmission, smoothed to a resolution R = 100, is shown for column densities 
N(R 2 ) = 10 16 , 10 18 , and 10 20 cm' 2 with the level populations of Fig. 1. 
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Fig. 4. — Fluorescent emission vE v per irradiated H2 molecule, for irradiation by I v oc v x l 2 with 
cutoff at 912 A (upper panel), or including Lyman continuum (lower panel). 



